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ABSTRACT 

Aims. We aim at measuring the stellar parameters of the two Chemically Peculiar components of the B9.5Vp HgMn + A2 Vm 
double-lined spectroscopic binary HD141556 (j( Lup), whose period is 15.25 days. 

Methods. We combined historical radial velocity measurements with new spatially resolved astrometric observations from PIO- 
NlER/VLTl to reconstruct the three-dimensional orbit of the binary, and thus obtained the individual masses. We fit the available 
photometric points together with the flux ratios provided by interferometry to consttain the individual sizes, which we compared to 
predictions from evolutionary models. 

Results. The individual masses of the components are Ma = 2.84 ± 0.12 M© and Mb = 1.94 ± 0.09 M©. The dynamical distance is 
compatible with the Hipparcos parallax. We find linear stellar radii of = 2.85 ± 0.15 Rq and R^, = 1.75 ± 0.18 Rq. This result 
validates a posteriori the flux ratio used in previous detailed abundance studies. Assuming coevaUty, we determine a sUghtly sub-solar 
initial metallicity Z = 0.012 ± 0.003 and an age of (2.8 ± 0.3) x 10^ years. Finally, our results imply that the primary rotates more 
slowly than its synchronous velocity, while the secondary is probably synchronous. We show that strong tidal coupling during the 
pre-main sequence evolution followed by a full decoupling at zero-age main sequence provides a plausible explanation for these very 
low rotation rates. 

Keywords. Stars: binaries: spectroscopic - Stars: cliemicaliy peculiar - Stars: individual: HD1 41 556 -Stars: rotation - Tecfiniques: 
interferometric 



1. Introduction 



> 



Chemically peculiar stars (CP), usually classified as Bp, Ap, 
HgMn, or Am stars, have challenged our understanding of stel- 
lar atmospheres for decades (Preston 1974). Although not com- 
pletely understood yet, the origin of the peculiariarities seems 
closely related to a slow rotation velocity (Abt et al. 1972) and 
to the presence of a magnetic field that could bring to the at- 
mosphere the stability needed for diffusion processes to be im- 
CN portant (Michaud 1970; Michaud et al. 1974). The CP stars of 
mercury-manganese (HgMn) type are exceptional cases where 
no clear magnetic field is detected but peculiarities (and some- 
I times spots) are still present on the stellar surface. Observations 
^ suggest a correlation between the presence of a weak magnetic 
field, abundance anomalies, and binary properties (SchoUer et al. 
r> 2010; Hubrig et al. 2012). More studies of HgMn stars are es- 
pecially important for understanding the role of the multiplicity 
and the magnetic field as well as the evolutionary aspects of the 
chemical peculiarity phenomenon. 

The abundance anomalies of the non-magnetic CP star 
HD141556 (x Lup) are among the most extreme of this class. 
This star provides a unique astrophysical light source for atomic 
spectroscopy because it is very bright and exhibits numerous 
very sharp lines in the UV and visible parts of the spectrum. It 
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has been observed at very high signal-to-noise ratio by the God- 
dard High Resolution Spectrograph on the Hubble Space Tele- 
scope and other optical instruments, with the goal to derive abun- 
dances that fully span the periodic table of elements. See for in- 
stance the paper series by Wahlgren et al. (1994), Leckrone et al. 
(1999), Brandt et al. (1999) Proflitt (2002), and Nielsen et al. 
(2005) about the so-called Lupi pathfinder project. 

Part of the complexity of the task is that HD141556 is in 

fact a double -lined spectroscopic binary (SB2) composed of a 
B9.5Vp HgMn primary and an A2 Vm secondary. In agrement 
with the properties of these classes, the primary seems not to host 
an organized magnetic field while the secondary does (Mathys & 
Hubrig 1995). On the other hand, this characteristic provides the 
chance to determine dynamically the masses of the individual 
stars and the distance to the system. With this is mind, this well- 
studied object becomes a perfect case to compare the stellar pa- 
rameters with the prediction from evolutionary tracks. However, 
absolute masses have not been determined so far, most proba- 
bly because of the rather small size (3 milliarcseconds) of the 
apparent orbit. This is the purpose of this paper. 

Section 2 presents the historically available radial velocity 
measurements as well as the new astrometric observations ob- 
tained with long-baseline interferometry. The data set is used to 
solve for the dynamical masses of each component and the dis- 
tance of the system. In Section 3 we constrain the linear radius 
of each component by fitting the integrated photometry and the 
new flux ratio in the H band obtained by interferometry. Results 
are compared with the prediction from evolutionary tracks, pro- 
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Table 1. Log of interferometric observations with PIONIER at the VLTI together with the measured astrometric positions of b (faintest) with 
respect to a (brightest) and flux ratio in the H band. The position angle (PA) is defined north (0 deg) to east (90 deg). 



Obs. JD 


Date 


Telescopes 


Spectral 


East 


North 


Separation 


PA 


Flux 


2400000+ 




configuration 


setup 


(mas) 


(mas) 


(mas) 


(deg) 


ratio 


55982.845 


2012-02-24 


DO-HO-Gl-Il 


3 channels 


+2.20 + 0.1 


+2.17 + 0.1 


3.09 + 0.1 


+45.4+ 1.8 


2.8+0.2 


55989.857 


2012-03-02 


Al-Gl-KO-Il 


7 channels 


-2.69 ±0.1 


-1.96 ±0.1 


3.33 + 0.1 


-126.0+ 1.7 


2.9 + 0.1 


55990.906 


2012-03-03 


Al-Gl-KO-Il 


3 channels 


-2.13 ±0.1 


-2.13 ±0.1 


3.01 ±0.1 


-135.0+1.9 


2.9 ± 0.2 


55992.875 


2012-03-05 


Al-Gl-KO-Il 


3 channels 


-0.15 ±0.1 


-1.44 ±0.1 


1.45 ±0.1 


-174.0 ±3.9 


3.0 ± 0.2 



viding an estimate of the age and initial metalUcity of the system. 

Section 4 discuses the state of orbit circularization and rotation 
synchronization of the system. The paper ends with brief con- 
cluding remarks. 

2. Masses and orbital determination 

Interferometric data were obtained with the PIONIER combiner 
(Le Bouquin et al. 201 1) and the four auxiliary telescopes of the 
Very Large Telescope Interferometer (VLTI, Haguenauer et al. 
2010). Data were dispersed over three or seven spectral chan- 
nels across the H band (1.50 - 1.80/im). Data were reduced 
and calibrated with the pndrs package (Le Bouquin et al. 201 1). 
Calibration stars were chosen in the JMMC Stellar Diameters 
Catalog' (JSDC, Lafrasse et al. 2010). A log of the observa- 
tions can be found in Table 1. The amount of observing time 
per night is about thirty minutes, calibrations included. The cal- 
ibrated square visibilities and closures phases are shown in Fig- 
ures A. 1 to A.4, together with the best-fit binary models. Thanks 
to the simultaneous combination of four telescopes, the binary is 
detected in all PIONIER observations with no ambiguity in the 
relative position of the two components. 

We used the LITpro^ software (Tallon-Bosc et al. 2008) to 
extract the binary parameters, namely the flux ratio (assumed 
to be constant across the H band) and the astrometric separation. 
The expected diameters of the individual components (» 0.4 mas 
and 0.25 mas) were taken as fixed parameters. This assump- 
tion has no significant impact on the results because these di- 
ameters are almost unresolved by the longest VLTI baselines. 
The uncertainties computed by LITpro are unrealistically small. 
We performed a bootstrap analysis by randomly selecting sub- 
sets of the observations and adding noise to these datasets. The 
uncertainty is given by the dispersion of the best-fit parameters 
when fitting these datasets. This procedure was repeated inde- 
pendently for each date. We obtained a typical dispersion of 
0.05 -0.1 mas on the astrometry. Consequently, we assigned a 
conservative uncertainty of +0.1 mas in both right ascension and 
declination. Results are summarized in Table 1 . 

Table 2 summarizes the radial velocity (RV) measurements 
extracted from Dworetsky (1972), which include older data from 
Campbell and Moore (1928). We arbitrarily assigned uncertain- 
ties to the RV measurements by looking at the dispersion of each 
individual data set in Figure I from Dworetsky (1972). 

The specificity of the present data set is that we have RV data 
for the two components of the binary and astrometric data of the 
relative orbital motion. A combined fit of these data allows one 
to determine unambiguously all masses and orbital parameters of 
the system. In a standard Keplerian model, the radial velocities 
i^rad,a and Vrad,b of the two Stellar components and the ix,y) rela- 

' cdsarc . u- strasbg . £r/viz-bin/Cat?II/398 
^ www.jmnic.fr/litpro 
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Table 2. RV data extracted from Dworetsky (1972) 



Obs. JD 


Vrad,a 


I'rad.b 


Error^ 


2400000+ 


(kms-') 


(kms-i) 


(kms"') 


18549.50'' 


+36.6 


-103.3 


5.0 


18556.52'' 


-59.5 


+52.0 


10.0 


19097.88'' 


+33.2 


-79.1 


5.0 


19099.87'' 


+37.7 


-100.1 


5.0 


19221.71'' 


+39.2 


-101.9 


5.0 


37458.809 


-51.51 


+43.79 


1.0 


37832.776 


+27.43 


-74.17 


1.0 


38429.071 


+38.34 


-94.19 


1.0 


38430.043 


+37.29 


-94.77 


1.0 


38447.022 


+17.18 


-64.57 


1.0 


38448.053 


-3.76 


-34.41 


3.0 


38804.035 


-70.25 


+59.17 


1.0 


33805.032 


-60.68 


+44.66 


1.5 


38838.977 


+10.92 


-59.39 


1.0 


38981.101 


+12.35 


-67.50 


1.5 



Notes. Uncertainties were assigned arbitrarily based on information 
available in Dworetsky (1972). Original RV data from Campbell and 
Moore (1928). 



Table 3. Best fit orbital elements and related physical parameters 



Parameter 


Value 


Unit 


Vo 


-16.06 + 0.26 


kms~' 




55.11+0.46 


kms"' 


Kb 


80.69 + 0.47 


kms ' 


h 


38434.4 + 1.7 


JD-2400000 


P 


15.256560 ± 0.000071 


days 


e 


0.0076 ± 0.0054 




i 


110.2 + 2.1 


deg 


n 


-127.1 + 1.5 


deg 




115.9 + 39.7 


deg 


^app 


3.259 + 0.067 


mas 


Reduced 


2.609 




^phy 


0.2030 ± 0.0029 


AU 


Ma 


2.84 + 0.12 


Mo 


Mb 


1.94 + 0.086 


Mo 


d 


62.27 ± 1.7 


pc 



five astrometric position (projected onto the sky) of component 



J.-B. Le Bouquin et al.: HD141556 



10 




-100 




3 4 
Mean Anomaly (rad) 





Aa(mas) 



Fig. 1. Best-fit orbit compared to the historical RV observations (left) and the new astrometric observations from PIONIER (right). 



b relative to component a read 

Vrad,a = + (cOS(w + v) + g COS Lj) , 

Vradb = Vq - K\, (cOS((y + v) + g COS (jS) , 



X — 

y = 
with 



Tapp (cos O cos((y + v) - COS / sin Q. sin(w + v)) 
Tapp (sin Q cos((y + v) + cos / cos Q. sin(w + v)) 



'app 



aapp(l - e ) 
I + e cos V 



(1) 
(2) 
(3) 
(4) 

(5) 



Here, K-^ and are the individual amplitudes of the radial ve- 
locity wobbles of both stars, vq is the systemic heliocentric offset 
velocity, e is the orbital eccentricity, and /, f2 and lo stand for the 
inclination, longitude of ascending node, and argument of peri- 
astron of the orbit with respect to a fixed referential frame. The 
OXYZ referential frame is chosen in such a way that the XOY 
plane corresponds to the plane of the sky, that the OZ axis points 
toward the Earth, and that Q. is measured counterclockwise from 
north (the OX axis is A6 and points toward north; the OY axis 
is Aa and points toward east), flapp is the apparent semi-major 
axis of the astrometric orbit in milli-arcsec, and v is the current 
true anomaly along the orbit, v is related to the time by standard 
Kepler formalism, via the introduction of the orbital period P (or 
the orbital angular velocity n) and the time of periastron passage 
tp. 

There are therefore ten independent parameters to fit, namely 
vo, K^, Kb, tp, n, e, i, Q., o), and flapp- Additional physical quan- 
tities can then be derived from the knowledge of these parame- 
ters. The RV amplitudes are for instance related to the individual 
masses Ma and Mb by 



K, 



a,b 



Ma,b n flphy sin ; 



Ma + Mb Vl -e^ 



(6) 



where flphy is the physical semi-major axis of the orbit (in AU). 
Thanks to our astrometric data, the inclination / can be fitted as 
an independent parameter, so that flphy can be simply deduced 
from A'a + K\,. Then, with the knowledge of n, the total mass 
Ma -I- Mb is derived from Kepler's third law, and the values of 
/Ta and A'b yield the individual masses Ma and Mb. Finally, 
the comparison between the physical (flphy) and apparent (flapp) 
semi-major axes gives access to the distance d of the system. 

We performed a combined Levenberg-Marquardt least- 
square fit of the data. The convergence toward a single and 



deep minimum is rapid and robust for a wide range of initial 
guesses, so that we are confident in our orbital determination. 
The error bars on the parameters were deduced from their co- 
variance matrix, which is relevant in the present case of a deep 
minimum where the errors are probably Gaussian. The result of 
this fit is summarized in Fig. 1 and Table 3. 

The very large time separation between the first RV obser- 
vations (1911) and the last astrometric measurements (2012) al- 
lows the orbital period to be determined with a very high ac- 
curacy (« 6 seconds, i.e., a; 5 x 10"^ relative accuracy). The 
physical separation corresponds to about 15 stellar radii of the 
primary. The orbit is compatible with a null eccentricity, be- 
cause the error bar on the eccentricity is comparable to the fitted 
value. Consequently, the w parameter is meaningless and uncon- 
strained. The estimated distance to the system is compatible with 
the Hipparcos parallax for this star 16.71 ± 0.27mas = 59.8 + 1 
pc (van Leeuwen 2007). The mass of the primary is slightly 
larger than the photometric estimate Ma = 2.72 from Guthrie 
(1986). 



3. Temperatures, sizes, and age 

We retrieved the observed spectral energy distribution (SED) of 
the sum of the two components using the tool VOSpec^ provided 
by the ESA Virtual Observatory (data coming from the follow- 
ing catalogs: lUE for the UV, Johnson and SDSS for the optical, 
2MASS for the near infrared, and IRAS for the infrared). In 
addition, the interferometric observations presented in the previ- 
ous section constrain the flux ratio to be 2.9 + 0.1 in the H band 
(1.63 jum). From spectroscopic analysis in the optical, Wahlgren 
et al. (1994) determined a slightly higher flux ratio of 3.92 at 
0.44 fim and 3.62 at 0.6 yum. All these measurements are re- 
ported in Fig. 2. 

We built synthetic SEDs for each component using the 
DFSYNTHE, KAPPA9, ATLAS9, and SYNTHE codes and tables (Ku- 
rucz 1991; Sbordone et al. 2004; Castelli 2005). We used the 
GNU/LINUX port from Sbordone"* and Castelli^ For the pri- 
mary, we furthermore computed a custom opacity density func- 
tion (ODE) table using the abundances published by Wahlgren 



www . sciops . esa . int/index . php?pro j ect=SATcSipage=vospec 
* atmos.obspm.fr 

' wwwuser . oat . ts . astro . it/castelli 
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Fig. 2. Observed SED of the integrated flux of the binary (gray points), 
modeled by a sum of two Kurucz models with the stellar parameters 
described in Sect. 3 (solid black line). The red and blue circles are 
the flux of the primary and the secondary obtained by combining the 
observed flux ratio with the total flux at nearby wavelengths. These 
points are lowered by half a decade for clarity purpose. The dashed 
lines represent the flux of the corresponding Kurucz models. 



stellar extinction in the solar neighborhood give an extinction 
of A(V) - 0.05, which has negligible impact on the results. 

The diameters and effective temperatures of both compo- 
nents were adjusted so that the predictions from the models 
match the integrated photometric measurement and the available 
flux ratio. We explored a grid of diameters from one to four solar 
radii and temperatures from 6 000 K to 15 000 K. The result in 
shown in Fig. 3. The confidence regions of our adjustment are 
represented by the solid lines. The results for the secondary are 
partially degenerated because the measurements of the flux ratio 
lie in the Rayleigh- Jeans domain. 

Several estimates of the effective temperatures and the sur- 
face gravities for this system exist in the literature, generally 
using Stromgren and Geneva photometry. Smith & Dworetsky 
(1993) measured the stellar parameters of the primary compo- 
nent and found = 10750 K and log(gJ = 4.0. Wahlgren 
et al. (1994) iteratively determined the parameters of the two 
components by a detailed spectroscopic modeling and obtained 
T, = 10650 K, Tb = 9 200 K, log(^,) = 3.89 and log(^b) = 
4.2. Hubrig et al. (1999) found = 10664 ± 400 K and 
log(ga) = 4.22 ± 0.05, when fixing the parameters of the sec- 
ondary to Tb = 9200 K and log(gb) - 4.2 (taken from Wahlgren 
et al. 1994). Our independent determination of the masses of 
each component allows us to convert the surface gravity g and 
the mass M into the physical radius R of the photosphere: 



R = 



GM 



(7) 



c 

DC 
g 

T3 

TO 

DC 




Fig. 3. Solid lines are the 3cr confidence regions (blue for the pri- 
mary and red for the secondary) when adjusting the photometry and the 
flux ratio with Kurucz models. The dot, dashed and dash-dot lines are 
the confidence region derived from the studies of Smith & Dworetsky 
(1993), Wahlgren et al. (1994), and Hubrig et al. (1999) respectively. 
The filled stars are the adopted stellar parameters. 



et al. (1994) before running flie ATLAS9 code. The surface grav- 
ity of the models were adjusted to correspond to the masses 
estimated in Sect. 2 although the impact on the result is very 
small. We used the distance estimated in the previous section 
to compute the apparent fluxes. Typical values for the inter- 



where G is the gravitational constant. The resulting confidence 
region in effective temperature and radii are represented by the 
dot, dashed, and dashed-dot lines in Fig. 3. When no other in- 
formation is given in the literature, we assumed an uncertainty 
of ±250 K on the temperature and +0.1 dex on the log(g). The 
error boxes overlay adequately with our confidence regions. 

To better constrain the stellar radii we decided to adopt 
the effective temperature of the primary determined from spec- 
troscopy by Wahlgren et al. (1994), i.e. = 10650 K. Doing 
so, we obtained a best fit with Tb = 9 100 K, R.^ = 2.85 ±0.15 Rq 
and /?b = 1-75 ± 0.18 Rq. Assuming an uncertainty of ±250 K 
on the eff'ective temperature of the primary does not significantly 
modify the stellar radii. The resulting modeled SED is displayed 
in Fig. 2, together with the available measurements. Our revi- 
sion of the stellar radii lead to a ratio /?a/^b = 1.63 ± 0.08. This 
is in reasonable agreement with the value Rn/Rh = 1.67 derived 
by Wahlgren et al. (1994) from the spectroscopic flux ratio. The 
flux ratio in the B, R, and H bands are all correctly reproduced. 
Since the latter is a direct and fully independent measurement, it 
validates a posteriori the flux ratio used in the abundance stud- 
ies from Leckrone et al. (1999), Brandt et al. (1999), and Nielsen 
et al. (2005). 

We then compared our estimates of the effective tempera- 
tures and stellar radii with predictions from evolutionary tracks 
computed for the masses determined in Sect. 2. This is shown 
in Fig. 4. The confidence boxes represent the intersection of 
the different constraints from Fig. 3. We used the tracks from 
Mowlavi et al. (2012) and a slightly sub-solar initial metallicity 
Z - 0.012. For both components, the confidence boxes overlap 
correctly with the evolutionary tracks. We explored how sensi- 
tive the agreement is to the chosen initial metallicity. We per- 
formed a grid of {Age, Z}, and compared the predicted radii and 
effective temperature to the observations. Doing so, we assumed 
that the two components are coeval and have the same initial 
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Fig. 4. Evolutionary tracks from Mowlavi et al. (2012) for slightly 
sub-solar initial metallicity (Z = 0.012). Dashed lines are isochrones 
for ages from 1 to 5 x 10^ years. The gray regions correspond to tracks 
compatible with the dynamical masses estimated in Sect. 2. The col- 
ored boxes are the effective temperatures and stellar radii discussed in 
Sect. 3. 

Table 4. Angular rotation velocities, assuming the radii determined in 
the present work and \\\nir values from the literature. 



Ref. 


Va sin ir 


Vb sin 








(kms-i) 


(kms-i) 


(jurad s"') 


(jurad s"') 


(1) 


1 


2 


0.5 


1.7 


(2) 


2.5 


2 


1.4 


1.7 


(3) 


2 


4 


1.1 


3.5 


Average 


1.8 ± 1 


2.7 + 2 


1.0 + 0.6 


2.4+ 1.8 



References. (1) Leckrone et al. (1999) ; (2) Hubrig et al. (1999) 
(3) Nielsen et al. (2005) 



metallicity. These are realistic assumptions for such a tight bi- 
nary of massive stars, whose probability to result from a capture 
is extremely low. The resulting map is shown in Fig. 5. Every 
initial metallicity in the range Z = 0.009 to 0.014 provides a cor- 
rect agreement between the tracks and the boxes, and respects 
coevality. The best agreement is formally found for an initial 
metallicity Z - 0.0115 and an age of 2.85 x 10** years. 



4. Circularization and synchronization 

The mechanisms responsible for the development of the chem- 
ical anomalies in the stellar atmosphere of HgMn stars are not 
yet fully understood. Slow rotational velocities are definitely the 
most important ingredients. For HD141556, the existing studies 
all found low small rotational velocities but disagree on the ex- 
act values (see Table 4). Measuring such a small broadening is 
indeed challenging. However, we note that the vsin/,- are still 
higher than the microturbulent velocities (^a = km s"' and 
= 2 km s"', Wahlgren et al. 1994). This makes us confident 
that the measured rotational velocities are, at least, reliable upper 
limits. The literature provides no uncertainties for these quanti- 



0.020 








0.2 


0.4 

1 


0.6 

1 


0.8 


1 



0.015 



N 



0.010 - 




Age (years) 

Fig. 5. Map of 'Logix^) versus initial metalicity (vertical axis) and age 
of the system (horizontal axis) when comparing the effective tempera- 
tures and radii to the prediction of the evolutionary tracks from Mowlavi 
et al. (2012). We used the dynamical masses estimated in Sect. 2 and 
assumed coevality and identical initial metallicity for the two compo- 
nents. 



ties and we used instead the dispersion of the existing measure- 
ments. 

It is worth questioning, then, whether the system has reached 
such low rotational velocities by means of synchronisation be- 
tween rotation and orbital motion, thanks to past or still present 
tidal interaction. To do this, we need to compare the true angular 
rotation velocities of the stars with the orbital angular velocity 
n - 4.7 yurad s"'. 

Tides in a binary system usually act to align and synchro- 
nize spins with orbital motion and at circularizing the orbit. It 
is well known that spin alignement and synchronization are usu- 
ally achieved before circularization. The virtually zero eccen- 
tricity of the pair is an indication for an already achieved tidal 
circularization, and thus spin alignement and synchronization. 
Consequently, we may reasonably assume that the stars have ro- 
tation axes perpendicular to their orbital plane (zero obliquities, 
= /). With this assumption and our knowledge of the stellar 
radii, it is then possible to derive the angular rotation velocities 
N-i and A^b of the stars: 



A^ = 



vsinir 
R sin / 



(8) 



The results are displayed in Table 4. An exact synchronism 
would imply N^-Nh- n. While the secondary could appear to 
be synchronized within the error bars, obviously the primary ro- 
tates sub-synchronously with respect to the orbital motion. Sub- 
synchronism is rare among short-period binaries, but some cases 
have already been reported (Casey et al. 1995; Fekel et al. 201 1). 
Strikingly, the occurence of sub-synchronism seems to be higher 
in binaries with an HgMn component (Guthrie 1986; Hubrig & 
Mathys 1995). Hence it is interesting to study the history of syn- 
chronization in the HD141556 system. 
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Table 5. Characteristic times for the different tidal mechanisms, com- 
puted at ZAMS. Details can be found in the appendix. 



Mechanism (ref .) 


^circ 




(^sync)^ 




(fsync)^ 




(yrs) 




(yrs) 




(yrs) 


Equilibrium tide (1) 


1.37 X 


IQi^ 


3.67 X 


10^3 


2.70 X 


Dynamical tide (2) 


5.86 X 


IQi^ 


1.23 X 


1012 


5.04 X 


Meridional circulation (3) 


1.94 X 


10i« 


3.20 X 


10^ 


1.79 X 



References. (1) Press et al. (1975) ; (2) Zahn (1977) ; (3) Tassoul & 
Tassoul (1992) 



There are different theories for tides in close binary stars 
that provide predictions for synchronization and circularization 
times. The first one is the equilibrium tide, where each compo- 
nent, distorted by the other, reaches an equiUbrium figure. Then 
the viscosity prevents the equilibrium bulges from permanently 
pointing toward each other. The result is a coupling between ro- 
tation and orbital motion. This theory was first described in de- 
tail by Alexander (1973). The synchronization and circulariza- 
tion times depend on the mean viscosity of the stars, which are 
poorly constrained parameters. But Press et al. (1975) showed 
that in radiative envelopes, turbulent viscosity dominates, and 
the authors deduced explicit formulas. 

Zahn (1975) introduced another mechanism called dynami- 
cal tides that mainly applies to stars with radiative envelopes. In 
these stars, non-adiabatic oscillations driven by one component 
onto the other are damped by radiative dissipation, which causes 
a torque that is responsible for coupling between rotation and 
orbital motion. Zahn (1977) gives characteristic times. 

A third mechanism was introduced by Tassoul & Tassoul 
(1992) and is based on tidally driven meridional currents within 
the stars. According to these authors, this purely hydrodynami- 
cal mechanism can naturally explain the high degree of synchro- 
nism and circularization observed among close early-type and 
late-type binaries. They provide characteristic times for synchro- 
nization and circularization that mainly depend on the masses, 
sizes, luminosity, and separation of the components. 

The characteristic times are summarized in Table 5, while the 
details of the computation can be found in the appendix. The dif- 
ficulty is always to estimate the various characteristic constants 
that appear in these formulas. Note that the actual circularization 
time is computed from the two separate circularization times 



^circ (^circ)a (^circ)b 

because the contributions of both stars to the circularization pro- 
cess are added. 

The mechanism proposed by Tassoul & Tassoul (1992) is by 
far the fastest, as already noted by Beust et al. (1997). However, 
the reality of this mechanism is questioned (Rieutord & Zahn 
1997). We also note that dynamical tides are more efficient than 
equilibrium tides. This agrees with Zahn (2008), who stressed 
that dynamical tides are the dominant mechanism in stars with 
an outer radiative envelope (which is the case here), and that 
equilibrium tides are stronger for stars with a convective enve- 
lope and for planets. 

What conclusion can we derive for the HD 141556 system? 
First, we note that irrespective of the mechanism, circularization 
is not expected to be achieved at the age of the system. This 
could appear to be surprising given the eccentricity we measure, 
but it is worth recalling that the strength of the tides is highly 



time-variable. Basically, for such stars, tides are expected to be 
much stronger in the pre-main sequence (PMS) phase, thanks to 
the presence of convective envelopes. Circularization is there- 
fore very likely to have enough time to be reached during the 
PMS phase. But in that case, the stars should also meanwhile 
have synchronized their rotations. Second, we stress that if the 
lOflssoul & Tassoul (1992) mechanism were active today, both 
lOftitrs would still be perfectly synchronized. The synchronization 
lOfmes associated to this mechanism are indeed very short. Hence 
any departure from synchronism should be quickly damped by 
meridional circulation. Thus the observed sub-synchronism of 
the primary (whatever its origin) indicates that meridional circu- 
lation is very likely inactive in this pair. 

Our proposed scenario is then the following: in the PMS 
phase, the tidal couphng of these stars was strong. It could be 
due to both dynamical and equilibrium tides. This strong cou- 
pling resulted in a rapid synchronization and circularization of 
the pair Then, after the zero-age main-sequence (ZAMS) and 
the settling of a radiative envelope, the tidal coupling dropped 
drastically. In the absence of any external perturbation to the 
orbit, the eccentricity remained low. But between ZAMS and 
today, both stars underwent a significant radius increase. With 
our knowledge of the current radii and of the radii at ZAMS 
(Fig. 4), we can compute the expected angular velocities for the 
present day A^now = n (/?zams/^now)^- This formula assumes per- 
fect synchronism at the ZAMS, conservation of the rotational 
angular momenta between ZAMS and now, and no change in 
the internal structures of the stars. We found A^a - 1.9/irads~' 
and A'b = 3.5 jurad s"' . These values are still slightly higher than 
that of Table 4, but, at first order, this simplistic scenario pro- 
vides a plausible explanation for the sub-synchronism of the two 
components. 

We stress that there is no definite theory for the explanation 
of sub-synchronous rotation. A classical scenario for slowing 
down stars is of course a magnetic field. The existence of an 
efficient magnetic braking has been proposed to explain why all 
the magnetic A stars (Ap) have very slow rotational velocities. In 
binary system, such a magnetic braking could overcome the tidal 
synchronisation and naturally lead to sub-synchronous rotation 
(Hubrig & Mathys 1995). This is especially true for HD141556 
because the tidal coupling seems weak according to the calcu- 
lation presented above. However, this scenario requires that the 
two components are, or were at some point, strongly magnetized. 
Moreover, it does not explain the almost perfect circularization 
of the orbit. 

5. Concluding remarks 

HD141556 is a promising target for a precision test of evolu- 
tionary models with chemically peculiar stars. Accurate radial 
velocity measurements (according to modern standards) and a 
few additional interferometric observations spread over the or- 
bital cycle should considerably reduce the mass uncertainties, 
which are currently at the » 5% level. A single interferometric 
observation in the visible could provide a direct estimate of the 
flux ratio in this band (existing estimates rely on spectroscopic 
analysis), which should help in constraining the temperatures 
and the ratio of radii. A direct measurement of the stellar di- 
ameters (0.43 mas) is within reach with the SUSl interferometer 
or a next-generation visible instrument at VLTl. 

Our proposed scenario to explain the current rotation veloc- 
ities impUes that the stars have been synchronized by tidal cou- 
pling in the final stage of contracting toward the zero-age main 
sequence, and then fully decoupled during their evolution on the 
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main sequence. Adding generality to this result, it supports the 
fact that many HgMn stars are found in young multiple systems 
and in young associations. It is conceivable that these stars have 
been synchronized and that their peculiarities were established 
as early as the pre-main sequence phase. A quantitative test of 
this scenario requires modeling the internal structure of the two 
stars at all times, including the pre-main sequence. 
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Appendix A: Interferometric observations 

Appendix B: Computation of chiaracteristic times 
for circularization and synchironization 

Appendix B.1: Equilibrium tide 

Press et al. (1975) provide explicit formulas. For each compo- 
nent of the binary, we have 
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Here A' is the angular rotation velocity of the star, g is its gyra- 
tion radius (the moment of inertia is defined as g^mR^), Rt is an 
effective Reynolds number that can be taken ~ 20, and K and 6 
are numerical parameters. According to Press et al. (1975), we 
can typically assume K = 0.025 and 6 -Q.l. The other symbols 
have their meaning introduced above. These formulas are valid 
for either component of the binary. In each case, mi stands for 
the mass of the star we are considering, and mi for the mass of 
the other star. 



Appendix B.2: Dynamical tides 
Zahn (1977) gave characteristic times. 
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The conventions are the same as above. E2 is a dimension- 
less constant characteristizing the strength of the dynamical tide. 
Zahn (1975) provided a table of E2 and gyration radius at ZAMS 

for various masses. We interpolated these values for the actual 
masses and computed the characteristic times, assuming n = N 
and e = 0. 



Appendix B.3: Tidally driven meridional currents 

According to Tassoul & Tassoul (1992), the characteristic times 
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L is the stellar luminosity and A' is a charateristic exponent 
that can be taken as 10 for stars with a convective envelope, and 
for stars with a radiative envelope, which is the case here. 
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Fig. A.l. Square visibilities (left), closure phases (middle) and corresponding reduced map (right) for the PIONIER observation of 2012-02-24. 
The corresponding best-fit binary model is shown in red. 
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Fig. A.2. Same as Fig. A. 1 for the PIONIER observation of 20 1 2-03-02. 
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Fig. A.3. Same as Fig. A. 1 for the PIONIER observation of 20 1 2-03-03 . 
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Fig. A.4. Same as Fig. A. 1 for the PIONIER observation of 20 1 2-03-05 . 
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